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We discuss the thermodynamic limit in the canonical area ensemble used in loop 
quantum gravity to model quantum black holes. The computation of the thermo- 
dynamic limit is the rigorous way to obtain a smooth entropy from the counting 
entropy given by a direct determination of the number of microstates compatible 
with macroscopic quantities (the energy in standard statistical mechanics or the 
area in the framework presented here). As we will show in specific examples the 
leading behavior of the smoothed entropy for large horizon areas is the same as the 
counting entropy but the subleading contributions differ. This is important because 
these corrections determine the concavity or convexity of the entropy as a function 
of the area. 
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I. INTRODUCTION 

The standard formalism of thermodynamics can be based on the study of the entropy 
as the fundamental function encoding the full thermal behavior of a system. This point 
of view was actually advocated by some of the founding fathers of the subject, in par- 
ticular by Gibbs [1]. The entropy is a function of the equilibrium states of the system. 
For example, in the case of a gas, the equilibrium states are described by the total energy 
E and the volume V of the container. Hence, the entropy is a function of (E, V). Other 
very important thermodynamical quantities are defined by derivatives of the entropy. In 
particular the inverse temperature and the pressure are given, respectively, by 

l_dS P _ 9S 
T ~ 9E ' T ~ 9V' 

The previous expressions only make sense if S(E, V) satisfies some obvious and necessary 
smoothness conditions. 

The entropy plays a central role in statistical mechanics. In fact, its determination by 
the counting of states subject to certain constraints (a fixed value of the energy, for exam- 
ple) can be taken as the starting point to understand the thermodynamical behavior of a 
system from the dynamics of its constituents. The actual result of such a counting usually 
suffers from a significant drawback: by its very nature it is always a staircase function of 
the natural variables describing the system (with jumps for the values of the basic vari- 
ables for which the number of allowed states changes). This means that the derivatives 
necessary to compute relevant quantities, such as the temperature, are either zero or fail 
to be defined. A well-known solution to this problem goes back to the classical paper by 
Griffiths [2] where the author shows how a smooth entropy -with the possible exception 
of isolated points that would mark the existence of phase transitions- can be found in a 
suitable thermodynamic limit (see [3] and also [4] for more details). This limit is defined by 
taking the size of the system to infinity while keeping some intensive parameters fixed. 
In the case of a gas, for example, the volume of its container is taken to infinity while 
keeping the number particle density and energy per particle fixed. In standard statistical 
mechanics texts this smoothing is sometimes carried out by considering discrete vari- 
ables as continuous (which amounts in practice to some sort of coarse graining). Though 
this is not a rigorously defined procedure in many cases it provides the correct answer. 

The identification of the laws of black hole mechanics [5] in classical general relativity 
suggested an unexpected connection between black hole physics and thermodynamics, 
in particular, a relationship between entropy and horizon area. The possibility of relat- 
ing the horizon area with entropy was proposed by Bekenstein [6] after the discovery by 
Hawking of the fact that the area of black hole never decreases [7]. The exact proportion- 
ality between area and entropy, S = A/4 (in units h = c = G=k B = l), was fixed after 
the discovery of the Hawking radiation [8] with a thermal spectrum. 

The study of the microscopic origin of black hole entropy, i.e. black hole statistical me- 
chanics, has received a lot of attention in the last fifteen years. The successful derivation 
of the Bekenstein-Hawking law S = A/4 has been claimed in several different approaches 
to quantum gravity and used by their proponents to highlight their relative merits. Two 
reasonably developed examples (but certainly not the only ones, see [9]) are provided 
by string theory and loop quantum gravity (LQG). In both cases there are schemes to 
define the entropy for certain types of black holes through the counting of microscopic 
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configurations. Both approaches boast successful derivations of the Bekenstein-Hawking 
law. 

The counting of states in string theory can be carried out for extremal [10] or quasi 
extremal black holes (see [11] for a recent review). Furthermore, the correspondence 
principle put forward by Horowitz and Polchinski [12] leads to a generic linear depen- 
dence between area and entropy for arbitrary black holes though the proportionality con- 
stant cannot be obtained by invoking only this principle. These computations rely on the 
asymptotic determination of the number of D-brane configurations in the limit of large 
charges obtained by using the Cardy formula. The resulting expression for the number of 
states can be compared with the area of a classical black hole carrying the same charges. 
In the relevant limit (i.e. for large charges) one indeed finds the sought for linear relation- 
ship between entropy an area with the 1 /4 proportionality factor. 

The loop quantum gravity framework approaches the same problem by modeling 
black holes with the help of the so called isolated horizons [13]. These can be understood 
as inner spacetime boundaries -with the topology § 2 x R- that satisfy some geometrical 
conditions related to the presence of marginally trapped surfaces and energy conditions. 
An important result in this setting is the possibility of deriving laws of isolated horizon 
mechanics that mimic some of the laws of thermodynamics [14]. The sector of general 
relativity consisting on space times with such inner boundaries admits a Hamiltonian 
formulation that can be taken as the starting point for its quantization [15, 16]. After in- 
troducing a suitable Hilbert space inspired in the ones customarily used in LQG (spanned 
by the so called spin network states) and solving the quantum constraint enforcing the 
isolated horizon conditions, it is possible to introduce a maximally degenerate density 
matrix p and compute the quantum entropy as Tr (p log p) . When this is done a linear re- 
lationship between area and entropy is found for large horizon areas. The proportionality 
coefficient is a function of the so called Immirzi parameter y that must be then suitably 
fixed to recover the 1 /4 coefficient of the Bekenstein-Hawking law [17, 18]. Although this 
is somewhat unsatisfactory, the choice of y is universal (in the sense that it is valid for 
all the types of black holed that have been studied) and hence the result can be used as 
a physical way to fix this otherwise undetermined parameter. It is important to point out 
here that the LQG formalism can be applied to physical (i.e. non-extremal) black holes. 

In the two settings considered above the counting entropy (referred to by some authors 
as statistical entropy) is a discontinuous function consisting in discrete steps. In order 
to get a suitable smooth function one has to consider the thermodynamic limit. It can 
be shown (by using, for example, a saddle point evaluation of certain integrals [19]) that 
the leading behavior of the true (smoothed) thermodynamical entropy coincides with the 
one corresponding to the statistical entropy. However, subdominant contributions can be 
different. This is specially important for black holes because the concavity or convexity 
of the entropy (related to the stability or lack thereof of the system) crucially depends on 
the behavior of these subdominant contributions. 

The purpose of this paper is to discuss the issue of the thermodynamic limit for black 
holes. As we will show this is not simply the large area limit of the counting entropy. 
A direct and concrete derivation of this limit in the string framework has not been per- 
formed despite the fact that a very detailed knowledge of the microstates for some types 
of black holes has been recently obtained (see [11] and references therein). The relevant 
counting of microscopic states in the LQG framework is also understood in great detail 
[20-22]. In particular the microcanonical area ensemble (the so called black hole degen- 
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eracy spectrum) is known in complete detail. Its Laplace transform defines the partition 
function in the canonical area ensemble; it can also be exactly obtained in this framework 
and used to illuminate some important features of the thermodynamic limit regarding, 
specifically subdominant corrections to the asymptotic value of the entropy as a func- 
tion of the area. In fact, this is the main goal of the paper. We will concentrate on the 
Domagala-Lewandowski (DL) [23] and the Engle, Noui, Perez (ENP) [18] proposals (see 
also [24, 25]) but our results can be extended to other LQG inspired models. 

The main conclusion of our analysis is that the subdominant corrections to the smooth 
entropy obtained in the thermodynamic limit differ from the ones corresponding to the 
counting entropy. This result is actually expected because the theorems that guarantee 
the existence of the thermodynamic limit for the area ensemble show that the entropy 
must be concave. However it is in apparent conflict with the asymptotic behavior of the 
counting entropy (that has subdominant corrections to the area law that are proportional 
to minus the logarithm of the area and are, hence, convex). 

The layout of the paper is the following. After this introduction we will briefly re- 
view in section II black hole and isolated horizon thermodynamics. We will then discuss 
in section III the introduction of the canonical area ensemble and give the correspond- 
ing partition function. We will show that, in the thermodynamic limit, the entropy is a 
smooth and concave function of the area and we will determine asymptotic behavior for 
large areas. As we will see the linear behavior is the same as for the statistical entropy but 
the logarithmic corrections change. An important comment that is relevant at this point is 
the fact that the interesting structure [26] found in the study of the black hole degeneracy 
spectrum cannot be present in the thermodynamical limit (though it may be relevant to 
study the detailed behavior of black holes in LQG and, in particular, Hawking radiation). 
We end the paper with our conclusions and comments and several appendices devoted 
to a brief review of the thermodynamic limit for some important sample systems and the 
asymptotic behavior of the partition functions. 



II. BLACK HOLE THERMODYNAMICS 

The classic no-hair theorems tell us that black holes are described by a very small set 
of physical parameters: the mass M, the angular momentum J and the electric charge Q 
(see, however, [27]). In particular the mass plays the role of the energy of a standard ther- 
modynamical system. For Kerr-Newman black holes we have that the horizon area A, 
the surface gravity k, the angular frequency of rotation Q. and the electrostatic potential 
O on the event horizon of the black hole are [28] 

A(M, J, Q) := ^ (j 2 + (M 2 + VM 4 - Q 2 M 2 - J 2 ) 2 ) , 
M(a/M 4 - Q 2 M 2 - J 2 ) 



k(M,J,Q) 
Q(M,J,Q) 
0>(M,J,Q) 



J 2 + (M 2 + v/M 4 - Q 2 M 2 - J 2 ) 2 
4ttJ 

M-A(M,J,Q)' 
47tQ(M 2 + v/M 4 - Q 2 M 2 - J 2 ) 
M-A(M,J,Q) * 
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Those quantities satisfy the first law of black hole mechanics 



dM = -^dA + OdJ + OdQ. (Ill) 

07T 



If we define the black hole entropy and the temperature by 

S(M, J, Q) := A(M ; QJ) = ^ (j 2 + (M 2 + V / ^W r F) 2 ) , 

TfM T Q) •= K(M >Q^ = 1 M(y/M^Q 2 M 2 -T 2 ) 

27t In J2 + (M 2 + v/M 4 - Q 2 M 2 - J 2 ) 2 ' 

the first law (II. 1) takes the usual form of the first law of thermodynamics 

dM = TdS + OdJ + OdQ , 

where 

1 9S O 3S <D 3S 



T 9M' T 3J ' T 9Q ' 

Notice that the terms Cld] and OdQ play the same role of — PdV for a gas, i.e. they are 
work terms. It is important to realize that the area does not play the role of the volume. 

Alternatively it is also possible to express the mass (energy) as a function of the en- 
tropy and describe the black hole states in terms of (S, J, Q). By using these variables and 
A(S, J, Q) = 4S we have now 



2\ V2 



M(S,J,Q) 



((Q 2 + S/7t) 2 + 4J 2 ) 
2{S/n) yi 

(SAr) 2 -(Q 4 + 4J 2 



K(S,J,Q)=27tT(S,J,Q)- 

2(SAr)^((Q 2 + S/7T) 2 + 4J 2 ) /2 

Q(S,J,Q) 



#(SJ,Q) 



(S/7r)V^((Q2 + S/7t) 2 + 4J 2 ) V2 

{Q 2 + S/n)Q 
(S/7r)V2((Q2 + s/7t) 2 + 4J 2 ) 



V2 



The temperature is defined, as usual, as the derivative of the energy with respect to the 
entropy at constant "volume" (i.e. constant J and Q) 

k _3M 3M _ dM 

Isolated horizons provide a generalization of the black hole event horizons. On one 
hand, the definition of isolated horizons is quasi-local, and it refers only to certain fields 
defined intrinsically on the horizon. On the other, only the intrinsic geometry of the 
isolated horizon is assumed to be time independent, whereas the geometry outside (not 
fixed by the one on the isolated horizon [29]) may be non-stationary. In practice this 
means that there are nontrivial examples of isolated horizons in addition to the event 
horizons of the globally stationary black holes (which are isolated horizons themselves). 
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As in the case of black holes, every state of an isolated horizon allows us to define the 
physical quantities (M, J, Q) or, equivalently, (S, J, Q). These numbers refer only to struc- 
tures intrinsically defined on the horizon, without any reference to the behavior at infin- 
ity 1 The first law of black hole mechanics can be generalized to the first law of isolated 
horizons and all formulas written above for Kerr-Newman black holes are still valid for 
(rigidly rotating) isolated horizons in the Einstein-Maxwell theory [13, 31]. 

It is widely accepted that a successful quantum theory of gravity must explain the 
previous results from a microscopic point of view, however, there are some difficulties 
(beyond the obvious one of the lack of a fully working quantum gravity theory). For 
example, if the a microscopic description of a large mass Schwarzschild black hole is ca- 
pable of explaining the relation S(M, 0, 0) = 47rM 2 , then the statistical mechanical density 
of energy microstates states has to satisfy log Q(M.) ~ 4nWl 2 for large M. This means that 
the canonical ensemble is ill defined because [28] 



Z(|3) = 



D(M)e- pM dM 



•oo 

2 



e 47 ™ V pM dM 



diverges for all values of (3. There are several possible ways to sidestep these problems. 
One can, for example, restrict oneself to working with the microcanonical (energy) en- 
semble. Another possibility is to put the back hole system inside a spherical cavity [32]. 
Finally one can follow the suggestions of [33, 34] and use an area canonical ensemble. 
This is the path that we will take here. 



III. BLACK HOLE AREA ENSEMBLES 

The purpose of this section is to define and discuss the microcanonical and canonical 
area ensembles in LQG. The idea of using these types of ensembles goes back to Krasnov 
[33-35] and is, somehow, a necessity in the LQG formalism as the counting of states is 
naturally done by using the horizon area instead of black hole mass (see [17, 36]). In the 
following we will restrict ourselves to the spherically symmetric case but our results can 
be extended to more general situations such as the ones discussed in [37]. A precursor of 
our work can be found in Meissner [38], where he solves the relevant counting problems 
by using Laplace transforms and, hence, he essentially derives the canonical ensemble 
(corresponding to the microcanonical area ensemble introduced by DL [23] for the com- 
putation of the statistical entropy according to the recipe given in [17]). Our presentation 
here ties a number of loose ends: 

• We show how the exact resolution of the combinatorial problems [20, 21, 39], that 
gives the microcanonical ensemble and the statistical entropy, is a direct conse- 
quence of the use of Laplace transforms. 

• We argue that the canonical ensemble gives a smoothed and well behaved entropy. 
This is necessary to have the possibility of using the standard formalism of thermo- 
dynamics. An important consequence of deriving the entropy from the canonical 

1 It is important to notice that, in contrast to black holes, the isolated horizon mass and angular momentum 
do not suffice to provide a characterization the time independent horizon geometry. This characterization 
is given in terms of an infinite set of multipoles [30] that capture the allowed distortions in the mass and 
angular momentum distribution on the horizon. 
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ensemble is the fact that the subdominant corrections to the entropy do not coincide 
with those corresponding to the statistical entropy. This result may be relevant out- 
side the realm of LQG and should equally apply to string theory inspired models. 

• The smoothed entropy is concave (i.e. the second derivative is negative) as a func- 
tion of the area. This is actually a very general result (see [2]) and is relevant to 
discuss the stability of the system. A word of caution may be necessary here be- 
cause standard black holes are unstable. The likely reason behind this discrepancy 
is the use of the area ensemble. Actually, we are not claiming that the black holes 
are stable, but rather show that the use of the area ensemble has important physical 
consequences. 

• The entropy vanishes for zero area. This suggest a version of the third principle in 
the case of black holes in LQG. 

• The entropy given by the area canonical ensemble for a single black hole corre- 
sponds to the thermodynamical limit for a ensemble of non-interacting black holes 
(similar in spirit to the Einstein crystal model discussed in Appendix A). 

The main approaches to the problem of counting the configurations giving the statis- 
tical entropy in the microcanonical area ensemble for a back hole are those of DL [23] 
and ENP [18]. For completeness we give the definitions of the statistical entropy in both 
cases. 

Microscopic black hole entropy: DL-counting. According to Quantum Geometry and 
the Ashtekar-Baez-Corichi-Krasnov framework, the counting entropy S^ icro (A) of a quan- 
tum horizon of classical area A, is given by 

S° L icro (A)=logO DL (A), 

where D DL ( A) is one plus the number of all the finite, arbitrarily long, sequences (mi , . . . , m N ) 
of non-zero half integers, such that the following equality and inequality are satisfied: 

N N ^ 

The extra one in the definition of I1 DL (A) comes from the trivial sequence. 
It is sometimes helpful to ignore the condition 



J^m-i = 



(the so called projection constraint) to get a simplified entropy S^ cro (A) = logO° L (A) that 
is useful to understand some features of the entropy in the LQG framework. This has 
the effect of changing the subdominant terms in the asymptotic behavior of the entropy 
and, hence, contains important physics. This prescription must be understood as a way 
to rephrase the original counting problem in [17] as one that can be solved by simpler 
methods. 
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Microscopic black hole entropy: ENP-counting (when y < y/3). The entropy S^ ro (A) 
of a quantum horizon of the classical area A is defined as 

sr cro (A)=io g a ENP (A), 

where n ENP (A) is one plus the number of all the finite, arbitrarily long, sequences 
(ji , . . . , )n) of non-zero half integers ji satisfying 

N 



2 



and counted with a multiplicity given by the dimension of the invariant subspace 

Inv((g)i[ji]). 

In principle the entropy defined above makes sense only for prequantized area values 
A = A£ s = 47ry£pk, k e N, however we extend the definition for arbitrary values of the 
area in the obvious way. In both cases the combinatorial problems involved in the com- 
putation of 0( A) can be exactly solved by using number theoretic methods and the result 
conveniently encoded in generating functions as shown in [21, 40-42]. In the following 
we will use units of 47ry£p 

In the thermodynamic limit, and irrespective of the model, the partition function per- 
particle Z for an ensemble of non-interacting "particles" can be computed as the Laplace 
transform of the corresponding number of microstates O of a single object (see Appendix 
A): 



Q(A) = Y_ D n 9(A - A„.) <£> Z(a) = a e" aA D(A) dA = Y_ D ne 



n=1 



,-aAn 



n=1 



where the integer numbers D n encode the black hole degeneracies associated with the 
area eigenvalues A n and 6 denotes the Heaviside function. Notice that the parameter <x 
is conjugate to the area and hence is not a temperature (which is conjugate to the energy). 
In the thermodynamic limit (see Appendix A), the average area is given by 

d 

a(a) = -— logZ(a) 
da 

and the (smoothed) entropy can be computed as 

&(<x) := <xa(<x) + logZ(a) . 

In practice, in order to express (and plot) the entropy as a function of the area, a i— > a(a), 
it is convenient to think of a as a parameter and consider the parametrized curve <x i— > 

(a(a),a(a)). 

In the DL scheme the partition functions Z DL (a) and Z° L (a) (with and without the 
projection constraint, respectively) can be read off directly from the integral expressions 
for the statistical entropy found in [21, 38]. When the projection constraint is not used 
Z° L (a) is given by 

Zf (a) = 1 -== . (III.1) 

1-2X~ =1 e-<Vw 
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If the projection constraint is incorporated we have, instead 

1 [ 2n dcu 



Z DL (<x) 



In 



1 - 2 L^ e^v^) cos cuk ' 



(III.2) 



As real functions of the variable <x both Z DL (cx) and Z° L (a) have singularities for the 
unique value (Xq l 6 IR satisfying 



1 -2^e- a oV*+ 2 ) = 0. 



k=1 



This can be directly seen in the case of Z° L (ct), whereas for Z DL (a) the integral in the 
auxiliary variable cu diverges if oc = c\q l = (0.746231 ■ ■ ■ ) (and hence Z DL (<x) is, itself, 
singular). This singularity controls the large area asymptotic behavior of the entropy 
whereas the asymptotic behavior of Z{oc) in the regime oc — > oo controls the limit a — > 
of the entropy. 



100 r 




FIG. 1. Plot of the smoothed black hole entropy in the DL case without the projection constraint. 
We also plot the dominant asymptotic behavior that reproduces the Bekenstein-Hawking law and 
the correction obtained by considering the first subleading terms log a+C*, with C* := 1 +log q° I i 1 . 
Notice that both the exact entropy and its asymptotic approximation, 0Lq L cl + log a, are convex. It 
is important to mention that no staircase structure appears and also that the entropy vanishes for 
zero area. The smoothed entropy is plotted by using the parametrization a i— > (a(a), &(a)) and 
numerically computing the partition function. 
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FIG. 2. Plot of the smoothed black hole entropy in the DL case when the projection constraint is 
taken into account (here C := \ + log q° L y 2 + \ log 2). The main features are the same as in Fig. 1. 
The leading behavior is the same but the coefficient of the logarithmic term changes. 

The asymptotic behaviors of Z° L (<x) and Z DL (<x) near <Xq l are, respectively, of the form 



a 

Zf(a)~-^ + £q-(a-cC) , \ « 



a 



a — oa 

n=0 



DL+ 




q 



DL OO 



Joe - cx? L ^— : 



where the (non-zero) coefficients 



n / 2 v a — a J , a — > a, 



DL+ 



n=0 



1 1 n DL 

n DL 1 n DL 1 / T* — 1 



can be obtained in a straightforward way as discussed in Appendix B. The asymptotic 
behaviors of the entropy for large values of the area are then given by 

of (a) ~ <x° L a + log a + 1 + log + 0(1 /a) , a -> oo , 

cr DL (a) - o# L a + - log a + - + log q° L 1/2 + ^ log2 + 0(1/a) a -> oo . 

By fixing the Immirzi parameter y as in [38] it is possible to recover the Bekenstein- 
Hawking area law. Notice, however, that there are logarithmic corrections that imply 
that the difference between the entropy and the Bekenstein-Hawking law increases with 
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FIG. 3. Plot of the smoothed black hole entropy in the ENP case. We have plotted the entropy 
only to areas of around 5 (in units of 47ty£p) because the linear regime is reached very early and 
the concavity would otherwise not be apparent. 

the area, i.e. lim a ^ 00 (ff(a) — oca) = oo. It is very important to notice that the corrections 
given above differ from the ones corresponding to the statistical entropy [21, 38] 

SSdcJA) ~ cCA + o(A), S- icro (A) ~ <A - \ log A + o(log A) . 

This difference is to be expected because the function A i— > <Xq l A — I log A is not concave. 

Plots of the entropy as a function of the area and the asymptotic approximations can be 
seen in Figs. 1 and 2. One can readily see that the entropy is a concave (negative second 
derivative) function of the area. It is also worthwhile to notice how fast the entropy 
approaches its asymptotic approximation (a concave function itself in the DL case). 

The partition function in the ENP case can be obtained from the results of [42], it is 
given by 

Z ENP fa) = - sin 2 cuduj 



71 



1 -2Xi e - a v /k ( k + 2 ) sin ( ( k + 1 ) to) / sin w 



The singularity <x^ NP = [0.86'\006 ■ ■ ■ ) in the partition function that dictates the large area 
behavior of the entropy is defined now by 



l-^Jk+l 



k=1 
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In the vicinity of this singularity the partition function behaves as (see Appendix B) 

oo 

and the asymptotic behavior of the entropy is 

a ENP (a) - < NP a + logqr + 0(1/a), a -> oo . 

As in the DL case, the vale of y can be fixed to recover the Bekenstein-Hawking [42]. 
Notice that for the ENP model there are no logarithmic corrections (only the constant 
term log ql NF ) and, hence, cr ENP (a) — a™ 15 a remains finite when a — > oo. 

In this case, the correction also differs from the ones corresponding to the statistical 
entropy [42] 

SZcro(A) - (Xq NP A - \ log A + o(log A) . 
IV. COMMENTS 

The main point that we make in the paper is that it is not enough to determine the 
statistical entropy of a black hole by counting its microscopical states but, in addition, 
one should go to the thermodynamic limit in an appropriate way. This is unavoidable 
in order to find a smooth entropy function that can be differentiated to compute phys- 
ical quantities according to the well established rules of thermodynamics. Smoothness 
is related to convexity properties that are relevant to understand stability in thermody- 
namical terms and define the single valued Legendre transforms used in the definition 
of thermodynamic potentials. It is only in this limit that the eventual singularities of the 
entropy can be used to identify and study phase transitions. These issues are relevant if 
one is going to be serious about black hole thermodynamics. 

An illustration of the kind of problems that can be expected has been discussed here 
in the context of entropy computations for black holes in loop quantum gravity. The 
discussion is facilitated by the availability of very detailed counting schemes that enable 
a precise counting of states as functions of the horizon area. Actually it is possible to 
obtain the (area) partition function for the system in closed form by using generating 
functions (as hinted in the paper by Meissner [38]). With the help of the partition function 
it is straightforward to derive the form of the entropy in the thermodynamic limit. 

The methods that we have used through the paper are based on very general results 
in statistical mechanics. This means that we can extend our work to other LQG inspired 
approaches in addition to the ones discussed in this paper. In particular, it is possible to 
do it for the modified flux area operator model with equally spaced area eigenvalues [43] 
with conclusions similar to the ones found here. The discussion presented in the paper 
should also be relevant for the proposals appearing in [44] and [45]. 

It is important to point out here that the thermodynamic limit is not just the limit of 
large areas. This is so for several reasons. 

• The large area limit as such does not lead to a well defined smooth or concave 
entropy function. A concrete smoothing procedure must be implemented. 
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• The role of the area as a thermodynamical variable for black holes must be carefully 
understood. In the fundamental example provided by the Schwarzschild black hole 
the correct way to interpret the results on the entropy is to state that it is a function 
of the mass (energy) S(M) = 47rM 2 . The horizon area is A(M) = ]6nM 2 and the 
temperature is T(M) = 1 /(87tM). It is true that the entropy is proportional to the 
horizon area (i.e. entropy and area are not independent) but, in order to get the 
temperature by the using the standard rules of thermodynamics the entropy must 
be expressed as a function of the mass. 

• It is not clear at all that the black hole area plays the same role as the volume of a 
gas. This means that it is not obvious what kind of intensive parameter should be 
kept fixed when defining the thermodynamic limit (number of punctures divided 
by horizon area?). This problem disappears (or at least takes a different disguise) 
if one considers an ensemble of "independent" black holes similar to the Einstein 
model for a crystal. Although this leads to a definite prescription to get smooth 
thermodynamical properties one would expect the actual microscopic gravitational 
degrees of freedom and their interactions to play a relevant role (as spin interac- 
tions do in a ferromagnet). Our examples are meant to illustrate the importance of 
considering the thermodynamic limit but the ignorance about the exact dynamics 
of the system prevents us from going further. 

• Finally as we have shown, the subleading corrections do not necessarily coincide 
for the smoothed entropy derived in the thermodynamic limit and for the statistical 
one. Care must be duly exercised then when comparing subleading contributions 
to the statistical entropy with their "macroscopic" counterparts (as is frequently 
done in the literature). 

We want to add several more comments. First we want to emphasize the fact that the 
mere use of density matrices to describe quantum black holes means that we are actually 
working with some kind of statistical ensemble. In practice we are forced to measure a 
large number of times on a system consisting of a single object and prepared according 
to some concrete prescription. Nobody would hesitate, in the context of quantum me- 
chanics, to talk about the density matrix for say a hydrogen atom as a way to encode the 
classical uncertainty involved in the incomplete preparation of such a system. The way 
entropy is defined in LQG relies on a "maximal degeneracy" density matrix and, hence, 
incorporates this type of uncertainty. The Einstein crystal-like model that we are using to 
introduce a thermodynamic limit is inspired by this point of view. 

There is a certain ambiguity in the definition of the number of microstates as a con- 
sequence of the alternative ways to think about the system. For example, in the case of 
LQG black holes one can wonder whether the different sets of punctures (with their la- 
bels) are the analogues of the energy levels of a harmonic oscillator or if they should be 
considered as the particles in a gas (in an analogous way, one can wonder if the different 
D-brane configurations for a black hole system in string theory are really its microscopic 
constituents or they should be thought of as the levels in a harmonic oscillator). 

As a final comment we want to mention the fact that one should really work with an 
energy operator instead of an area operator. Although this is not available in the LQG 
models that we are using, one could try to postulate it [employing, for example the stan- 
dard one for Schwarzschild black holes M = y/X/ {4yfn)] and deal with the lack of a well 
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defined canonical ensemble by working directly with the microcanonical one in order to 
go to the thermodynamic limit. The results derived with the canonical area ensemble 
cannot be directly generalized to get the solution to this problem but on the other hand 
we expect that the general results that we have found (in particular the different sub lead- 
ing behaviors of the statistical entropy and the smooth entropy in the thermodynamic 
limit) will be generally true. 
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Appendix A: The thermodynamical limit for simple systems. 

In this appendix we briefly review the thermodynamic limit for some important sam- 
ple systems (see [2-4, 46] for more details). 

The quantum description of a N -particle system confined in the domain A C t 3 is 
provided by a Hamiltonian operator H(N; A). Let 



Q(E,N;A) =Tr(e(E-H(N;A))) 



be the number of eigenstates of H(N; A) not exceeding E. The counting entropy in the 
quantum microcanonical ensemble is defined as 

S micro (E,N;A)=logQ(E,N;A). 

As a function of the energy (taken as a real variable) it is a staircase function. This means, 
in particular, that it is either non-differentiable or the derivative is zero, a fact that pre- 
cludes the use of the standard formulas in thermodynamics. In the thermodynamic limit, 
the entropy per-volume is given by 

S micr0 (eV(A),pV(A);A) 
a(e, p)= lim — — , 

V(A)->oo V(A) 

where V( A) is the volume of the domain A, e is the energy per volume and p the number 
density. In some cases, the role of V( A) can be played by other extensive parameters (for 
example by the number of particles itself). With the exception of a set of zero measure 
(that corresponds to phase transitions) the function cr satisfies some regularity properties, 
in particular it is concave as a function of the energy e and, consequently, it is differen- 
tiable. This allows us to define state functions such as the (inverse) temperature 

_., , 1 9o\ , 
T[e,pJ 3e 



In the quantum canonical ensemble the temperature characterizes a thermal bath with 
which the thermodynamical system interacts (exchanging energy). In this scheme the 
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fundamental object that contains the relevant information about the system is the parti- 
tion function 

Z N (|3,A)=Tr(exp(-|3H(N;A))J 
or, equivalently, the free energy 

F N (|3,A) =logZ N (|3,A). 

In the thermodynamic limit, the free energy per volume is 

r( n , r W)(P,A) r logZ N(A] ((3,A) N(A) 
f(p > p)= i™ V (A) = V(A) ' Wh6re 9 = V(A) * 

In this scheme the energy per volume is a derived quantity given by 

3f 

e(P,p) = -gp(P,p) 

and the entropy per particle is obtained by a Legendre transform 

a((3,p) = (3e((3,p) + f((3,p). 

The existence of the Legendre transform is guaranteed now by the convexity properties 
of the free energy defined in the thermodynamic limit. Finally it is important to point 
out that for reasonable interactions the thermodynamic limit is the same for both the 
microcanonical and canonical ensembles, in particular we have 

a(e,p) = a(|3(e,p),p) . 

We give now some examples relevant for the discussion in the main body of the paper. 

Example 0.1 One-dimensional Einstein crystal. 

Let us consider a system consisting in N non-interacting one dimensional oscillators. The 
Hamiltonian is then H(N) = ©^H,, where each Hi = am*a.t is a (normal ordered) one- 
dimensional harmonic oscillator. For this system, it is straightforward tho show that in 
the quantum microcanonical ensemble 

^L(vrM N+ r). 

n=0 v 7 v 7 

S m icro(F, N) = logO(E, N) = log ( N + J^ E/a,J 

It is important to notice that, irrespective of the number of oscillators N e N, the counting 
entropy E i— > S micro (E, N) is a staircase function of E e [0, oo). However, in the thermody- 
namic limit the thermodynamic entropy e i— > cr(e), defined as 

S micro (N £> N) 1 /N + LNe/cuJ 

a e = km = lim — log 

= ^log(l + cu/e) + log(1 + e/aj) , (A.1) 
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is a smooth, concave, function for e G [0, oo) that satisfies a(0) = . The temperature is a 
derived quantity 

da logO+cu/e^ 
de cu 

On the other hand, if we take the thermodynamic limit within the canonical ensemble 
point view, the thermodynamic free energy per particle f coincides with the logarithm of 
the partition function of a single quantum harmonic oscillator, i.e. 

m = Hr = = lo S (Tr(exp(-PH,))) = log (f_ = log (j^) 



In this approach the mean energy per particle, e( (3), is 

df , „ , cu 



d|3 ' eP"- 1" 
Finally the entropy is computed as a Legendre transform 

&(P) = (3e(|3) + f ((3) = - log(1 - e^) . (A.2) 

As expected, (A.l) and (A.2) are equivalent 

a(e) = ff(P(e)), ff(P) = a(e(P)). 

An interesting fact, relevant in this example and due to the non-interacting nature of its 
constituents, is the possibility of finding the explicit relationship between the true and the 
statistical entropy for a single object. If we do not take the thermodynamic limit in the 
microcanonical ensemble and consider the N = 1 case, denoting Oi(E) := D(E,N = 1) 
and S^ (E) = S micro (E, N = 1 ), we get 

Oi(E) = expSi(E) = ( ] + LE/a, -h =1 + LE/cuJ = f^9(E/cu-u) 

^ ' n=0 

The Laplace transform L of exp Si, as a function of (3, is given by 

1 00 11 

£(Cl h p) = £(ex pSl , p) = - Y_ ^ = p • 

n=0 

Hence, as expected, the canonical partition function of a single quantum harmonic oscil- 
lator is, essentially, the Laplace transform of the exponential of the counting entropy in 
the microcanonical esemble 

Zt(P) = P£(a b P) = p£(expS b P) . (A3) 

The free energy per particle in the thermodynamic limit in the canonical ensemble is 
given by 

f(P)=logZ 1 (P)=log(P£(expS 1 ,P)) 
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and 

e(|3) = --^log(|3£(expS 1> |3)) 
Finally the relationship between both entropies is 

&(P) = -P^logO ( C(expS 1 ,P))+log(P£(expS 1 ,P)) 
Example 0.2 One-dimensional periodic lattice (phonons). 

As it is well known, (see, for example, [46]), the normal modes of a one-dimensional 
periodic lattice are described by the (normal ordered) Hamiltonian 

N 

H(N)=^a) n (N)<a n 

n=1 

where 



tu n (N) = 2a> 



sin 



and c is the lattice constant. The partition function is given by 



z N(p)=n 



i 



n LAl-exp(-(3cu n (N)) 
and hence, in the thermodynamic limit, the free energy (per unit length) is 

1 



1 1X1 / 
f((3)= lim Vlog — 

N^oo NC VI — 

n — 1 >■ 



n=1 

1 



exp(-|3o) n (N)) 



In 



r 2n/c / -J 

o l0g Vj -exp(-2|3o> sin(ck/2))J dk " 



Example 0.3 Black body radiation (photons) 

In this case, the extensive parameter that we will use in the thermodynamic limit is V = 
L 3 . The Hamiltonian is now 



H(L) = ^^a) n (L)a; A (L)a n ,A(L), 



■ A=± 



where A — ± are the two polarizations of the photon, n = (ni , n 2 , n 2 ) G Z 3 and 



2n 



uoM := VWV k n (L) = T n 



The partition function for finite volume is the product of the partition function corre- 
sponding to each normal mode 



* (W -n( i--»(-*^) )' 
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In the thermodynamic limit the free energy per volume is 

2 «- . / 1 



L— >oo 



m = Hm = Hm - Y log - 

L^oo L- 5 * \ 1 



liny 



V 
log 



nel 3 



1 



1 _ e -pVj? 



d 3 k = 



exp ( - (3oj n (L)) 

poo 

log 



1 



7t z 



1 



1 - e-$ w 



cv 2 da) = 



71 



45(3 3 ' 



and the energy density (energy per unit of volume) of the radiation is given by 



df, n , 1 f°° tu 3 dtu 



71 



o e^-1 15(3 4 



The entropy as a function of the temperature is 



&((3) = |3e((3)+f(|3) = 



47t 2 



45(3 



3 ' 



hence, the entropy is a concave function of the energy 

4 7T 1 / 2 



o-(e) = 



,3/4 



Appendix B: Asymptotic behavior of the partition functions 

Let us show first that 

DL 00 



n=0 



As Z° L (oc), defined in (III.l), is a meromorphic function the relevant coefficients can 
be simply obtained by computing limits when a — > ajf. In particular, to the order that 
we are considering in the text of the paper, we only need to know 



DL 



lim fa-cc° L )Z° L fa) = 



1 



The expression for Z DL (<x) given in (III.2) involves an integral so, in this case, we cannot 
proceed as before and, in fact, we actually need to find the asymptotic behavior of the 
integral in the limit a — > a° L+ . The best strategy in this case consists in writing 



Z DL (ct) 



1 

2~7T 



doj 



with 



oo oo oo 

Q DL (cc,o>) := l-2^e- a v / ^+ I 'cosa)k= (l -2 Y_ e^v 7 ^ 2 ') + (4 Y_ * 

k=1 k=l k=1 



sin — ) 
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We define now 



oo oo 

Q° L (cc, cv):=(]-2Y_ e - a v / ^) j + (4 kV^*^ s in 2 

k=1 k=l 



CU 



k=1 

-2r„\ , T32i 



C 2 (a) + B z (a) sin 



CU 



and 



G DL (cc, cu) := Qg L (a, cu) - Q DL (a, cu) = 4^ e^v 7 ^) (V sin 2 ^ _ sin 2 ^ ? 

k=2 



to write 



\DL 



+ 



1 



1 G DL 

+ 



This leads to the following expansion 



1 



Q8V Qo L Q DL Q 



DL 




Q DL Q^feVQ? 

This expansion is uniformly convergent if oc > ocq > o$ l and for all cu e [—71, 7t] because 
|G dl /Qq l | < Co < 1 (where Co ~ 0.87). As a consequence of this we can write 



Z DL (cx) 



1 

2?t 



dcu 



n Q^(a,cu) + 27t^J 



71 (G DL (oc,cu)) k 



dcu . 



(B.l) 



Only the first integral diverges when a — > <Xo L+ and, in fact, it can be computed in close 
form by using 



1 ; 

2~7t 



dcu 



1 



- n C 2 (cc) + B 2 (cx)sin 2 (cu/2) |C(ct)| y/C 2 {oL) + B 2 {a) * 



When a — > a? L+ we have that 



C(a) = 



and, hence, 
1 rn 

27t 



dcu 



+ ^c n (a-cx D T +V2 , 



n=1 



n C 2 (ct) + B 2 (a)sin 2 (cu/2) |B«-)| (a - ag 1 



n=0 



In the previous two expressions c n and d n are real coefficients. The remaining integrals 
in (B.l) can be seen to have the form 



1 1 
In 



1 



(Qo^cu))^ 



(sin 4k Y. Snflc, a)Pn( sin 2 y) dcu (B.2) 



20 



for some regular functions g n and polynomials P n . This can be seen by expanding 
~ 2 ^T = ^ -coskcu) = 1(1 -T k (cosa))) = l(l -2T k (l -sin 2 ^)) 



sin 



in terms of the Tchebycheff polynomials T k and using the fact 2 that the lowest degree 
monomial of (1 — T k (1 — 2x 2 ))/2 is k 2 x 2 . We need to use now 



J_ 
In 



sin 2n (a)/2) 



-duo (B.3) 



(C 2 + B 2 sin 2 (cu/2))™ 



k=0 



B 2 ^^VV ~ V 2C 2 y V 2 ' 2 ' '(2C 2 + B 2 ) 2 



Notice that, in our case we will always have n > m owing to the presence of the term 
sin 4k (cu/2) in (B.2). The hypergeometric functions appearing in the previous expression 
are of the form 

F N (z) := 2 Fi —z— , l;z 



2 ' 2 

It is straightforward to see that F N is a polynomial in the variable z for N < and satisfies 

Fn(z)= (1^-V 2 - for N = 1 » 2 ,... 

It is possible to show that (B.3) are always analytic functions in C in a neighborhood of 
C = when n > m. Indeed, if the index k in the sum satisfies k > m the corresponding 
term in the r.h.s. of (B.3) is the product of C' 2n k ' and an analytic function in C because the 
hypergeometric function is a polynomial. For the terms with m > k the hypergeometric 
function behaves as 1 / C 2(m_k)_1 and hence the full summand is of the form C 2 ( n - m ) +1 times 
an analytic function in C. We then conclude that (B.3) are analytic functions of a/oc — a° L 
and then Z DL (a) must have the form 



,,DL OO 

^Zlll L V" r, DL (n, — ^DLW 2 



with 



DL+ 



9-1/2 _ 



|B« L )| 

In practice some of the next order terms in the expansion for Z DL (a) can be computed 
directly without using the series expansion introduced above. For example it is straight- 
forward to see that 



n DL — — 



' dcu . 



2 This can be proved by using the following explicit expansion 

Ln/2j 



T k (x)= Y_ f^W-l)**' 

k=0 ^ ' 



k^n-2k 
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The asymptotic behavior of 

1 rn 



Z ENP (cx) = 



71 



sin cu dcu 



- n 1 -Iw e- tx v / ^+2) s in((k+ l)cu)/sina) 



can be discussed along similar lines so we give only the most relevant steps. In this case 
we can write 



Z ENP (<x) = 



1 rn 



7T 



sin a> 



-dcu 



with 



k=1 



Q ENP (a, cu) 

Q E -(a,a,) := 1 - f Sm(U1)( "e-^. 
tt sincu 

We introduce now 

/ oo 

Qr(a,cu) := l-^(k+l)e- 



- av /k(k+2) j 



k=1 



+ -^k(k+1)(k + 2)e 



-«VkO^) \ sin 2 



CU 



k=1 



and 



G ENP (cx,cu) := Qncx,o))-Q ENP (a,a>), 

which, as in the case of Z DL , is actually proportional to sin 4 (cu/2). Now we can write the 
expansion 

1 ^ (G-(«,u,))W a , 



and show that it is an analytic function of a/ oc — af F by following the same procedure 
that we used for Z DL . Though the previous series provides a way to compute the coeffi- 
cients of the power series for Z ENP (<x) it is better, in practice, to compute them directly In 
this case we have, for example, 



ENP _ 

Mo — 



1 



7T 



sin cu 



~ n YZU ((k + 1 ) sin cu - sin(k + 1 )cu) e-^v 7 ^) 



dcu 
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